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ABSTRACT 


Field-aligned electric currents flow in the magnetosphere in many situations 
of fundamental geophysical interest. It is shown here that the incoherent back- 
scatter technique can be used to measure these currents when the plasma line 
can be observed. To our knowledge this has not been proposed before. The 
technique provides a ground-based means of measuring these currents which 
complements the rocket and satellite ones. 
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INTRODUCTION 


Field-aligned o»’ ents in the magnetosphere exist in various geophysical 
contexts of fundamental interest. In the setting of auroral physics they were 
originally postulated by Birkeland and are now widely known by his name. 

From the experim ntal viewpoint most is known about these current in the 
auroral setting from the work of Zmuda et al (1966, 1967, 1970) Choy et al (1971), 
Cloutier et al (1973, Fairfield (1973), Thiele and Praetorious (1973, Vondrak et 
al (1971), Kawasaki et al (1974) Zmuda and Armstrong (1974), Ledley and Farth- 
ing (1974) and others. From the theoretical point of view these auroral currents 
have appeared in the works of Alfven (1939, 1940), Piddington (1959), Cole (1961, 
1963, 1971, 1974), Fejer (1963), Kern (1962), Bostrom (1964), Cummings and 
Dessler (1967), Swift (1963, 1965), Shield et al (1969), and others (see Anderson 
and Vondrak 1975, for a review). In the context of the quiet variations of the 
geomagnetic field they have appeared, e.g., in the work of Dougherty (1963), Van 
Sabben (1962), Mishin et al. (1968), Lemaire and Scherer (1974), Cocks and Price 
(1969), and Stening (1969). The inference of field-aligned current strengths using 
magnetometers on board satellites is made difficult by choice of baselines in the 
records and also by the complications of temporal as well as the spatial varia- 
tions of the magnetic field. 

As will be demonstrated the incoherent backscatter technique is capable of 
giving absolute measures of the field-aligned current from the strongest known, 
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to ones which are two orders of magnitude lesr,. This constitutes a great im- 
provement in our potential for knowledge of these currents. 

The field aligned current density J is given by 

ill s Jin + J.n = n 8 e(v t -yJi, W 

where 

v. = velocity of ions, electrons 
1 , © 

e = electronic charge 
n 0 = electron density 

II denotes the component along the magnetic field (B) 
Consider an incoherent s ockscatter device receivi"' 1 ' signals from the direction 
of the magnetic field. It is well known (see Bauer 1975 for a review) that tt e 
can be found from the intensity of the signals received; v i(J can be found from 
the shift of the so called "ion spectrum"; and we show here that v 0 |, can be 
found from the shift of the "plasma line" provided it is excited, e.g., by photo- 
electrons or auroral electrons at the same time. It follows that the current j| ( 
can be measured absolutely and directly with this ground-based technique. 

THEORY 

Measurement of v f| | 

Since the theory of measurement of n e and v l( | is well documented we give 
here the theory for the measurement of v e]t . Since we shall be concerned only 
with velocities h to B we shall hence forth drop the subscript || . It is supposed 
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that photoelectrons or auroral electrons traveling along the geomagnetic Held 
excite longitudinal waves in the ionospheric plasma. 


Consider then the dispersion relation for longitudinal waves of angular fre- 
quency (w) in a plasma without magnetic field, or for a wave vector k parallel 
to the magnetic field. It can be written (Stix 1962) 



f 0 (v) d v « 

+ bp 

(a) - kv) 2 P 


fj(v) dv 
(b) - kv ) 2 


( 2 ) 


where f B (v) is the one dimensional electron velocity distribution, f A (v) the one 
dimensional ion velocity distribution, <y pe the electron plasma angular frequency, 
<y pl the ion plasma angular frequency. Note that f e (v) relates to both the thermal 
and the suprathermal electrons. 


In the present problem we are only interested in thj real part of the solu- 
tions in the vicinity of ± a) . In this case equation (2) can be approximated by 
a power series expansion: 


<’?>' - 3) 


In order to simplify this expression it is necessary to make order of mag- 
nitude computations for co in the vicinity of ±w pe , The dominant term of the 
ionic component is, for oxygen ions, approximately equal to 3 X 10“ 5 times the 
major electronic term. The term dependent upon the mean square of the elec- 
tron velocity is of the order of 1/10 of the electron dominant term for mean 


3 


Ionospheric conditions (T e 1000K, f p 5 x 10 6 Hz) and a wavelength of 0.3 m cor- 
responding to a working frequency of 10 9 Hz, The term dependent upon the mean 
electron velocity roaches 1/100 of the dominant term for the same ionospheric 
conditions and a mean velocity of 5000 m/s. 


Therefore the ionic component can be neglected with respect to all the elec- 
tronic components and the dispersion relation for the solutions in the vicinity of 
rfc&j can be written as: 

pe 


or 


1 = 


In the vicinity of 


(a> - k <v o )0 2 
± o) p0 this expression has two solutions: 

" k < V e> = ± W 0 
"r = k < V o> ±a >0 


1 + 


3k 2 


[ k < v 0 >) 


O'. -<Ve>) 


2 


(4) 


(5) 


In other words the two solutions are now located symmetrically with respect 
to k<^ ( rather than with respect to zero. These solutions are approximately 
given by 

t(3KT,)k 2 (6) 


where T e is the apparent temperature of the electron gas (including the photo- 
electrons), and K is Boltzmann's constant. 

Whatever is the second term of the right hand side of (6) the important point 
is that a drift of the whole electron gas induces a doppler shift in the location of 
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the two solutions of the dispersion relation close to ±w pc , in terms of an inco- 
herent scatter experiment a drift of the electron gas Introduces a summed shift 
in the locations of the two plasma lines (1, 2) with respect to the central (i.o. 
transmitted) frequency (w c ) 

= 2k O 

e.g. for an electron drift of 1000 m/s, and a wavelength of .30 m 
A(J« r - = 2 x x 2 x 10 3 rad s" 1 


« 1.3 x 10 4 x 2tt rad s -1 

which can be measured. It follows that 

Am. _4V„ _4J e|| 

cj eo c n ec 

c c e 



The measurement of Aw /w, is then convertible to a measure of v 0 . 

Errors 

As for the ion spectrum the main cause of error is of statistical origin. 
However, contrary to the case of the ion spectrum the spectral shape of the 
plasma line is not determined by the ionospheric parameters alone but also by 
the distribution of the electron concentration in the volume being observed, i.e. 
by the geometry of the experiment. We make the reasonable assumption, in view 
of current experiments, that the plasma lines can be sketched as shown on Fig- 
ure 1. The flanks of the lines are characterized by approximately linear slope 
from the maximum signal S mQX to the noise level over frequency intervals A f. 
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We assume further that filters of width Af are being used. This leads to the 
following estimate of the error ca the electron velc jity for a noise level N 
larger than the signal (Petit 1968, Evans 19G9) 


Av o 


N v'Af' k 

s miix /r /2 


m/s 


where t is the time of integration k the wavelength. Table X gives the errors 
in velocity v e or electron current density (for an assumed electron density of 
5 x 10 ll m" 3 ), for A f varying between 25 and 100 kHz, a signal to noise ratio 
varying between 1/100 and 1, an integration time of 1000 sec and 'k ~ 0.3m. It 
will be shown that the accuracies mentioned in Table I are suitable for most cases 
for the measurements of field- aligned currents during quiet or disturbed conditions. 

In Figure 2 are shown values of Aw /w c as functions of j o( | and n 0 . With 
presently available equipment, values of Aw /w, 0 f order 10~ 6 are easily ob- 
tainable. On this figure are indicated ranges of current densities appropriate to 
disturbed and quiet conditions. Also indicated on Figure 2 are values of v 0 for 
various values of j e n and , As seen on Table I signal to noise ratio of 0.01 
should allow measurements of currents (j e || ) for disturbed conditions, and signal 
to noise ratio of 0.1 should allow measurements of j 0 || in both quiet and disturbed 
conditions. 


The ion drift velocity v. y is obtained from the shift of the "ion spectrum." 
The two measurements combined according to equation 1 will produce the field- 
aligned current j u . 
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Several special oases suggest themselves (a) ^ = 0, combined with bulk 
flow of plasma along B. In this case both the "ion spectrum" and the "plasma 
lines" are doppler shijf&jd by equal amounts. 

(b) j o || ^ 0 but jjn = 0. In this oase the "ion spectrum" will be unshifted 
but unsymmetrioal on account of the electron motion, and the plasma lines will 
be doppler shifted. 

(c) In the general case j o( | 4 0 and ^ 0. The ion-spectrum will be- 
shifted and unsymmetrioal and the plasma lines doppler shifted. The asymmetry 
of the ion spectrum can be related to - v o ). However, the asymmetry of the 
ion spectrum would relate only to the relative drift of the thermal electrons with 
respect to the ions and might indicate, if ever measurable, a relative drift of the 
photo electrons or auroral electrons with respect to the thermal electrons by 
comparison with the shift of the plasma lines. This will be particularly important 
in deciding the contribution of energetic and thermal electrons to Birkeland cur- 
rents in the auroral zone. In situations where the photoelectron contribution to 
the current is negligible compared to that of thermal electrons, this will provide 
a redundant check on the electron velocity measured by the shift in plasma lines. 

GEOPHYSICAL SITUATIONS OF INTEREST 

Electric currents flowing along the geomagnetic field occur in many geo- 
physical situations of interest. However, it is difficult to measure them in situ. 

Of the ground based techniques, incoherent scatter offers the opportunity to study 
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these currents in considerable detail both in tlxno, and altitude distribution. Wo 
present a disoussion of the various situations in which the currents occur. 

In auroras 

Birkeland considered field-aligned current in relation to the auroral 
eleotrojet. Theories have been proposed for them in this context. Measure- 
ments from satellite borne magnetometers have been interpreted by Zmuda and 
Armstrong (1974) in terms of field aligned currents of strengths up to about 
5 x 10“ s A/m 2 . Such a current requires a charge flux of 3 X 10 14 particles/m 2 S. 
With ionospheric densities of about 5 x lo 11 /m 3 , this implies (in the case where 
all the current is carried by electrons) an average electron velocity of 6 xl0 2 m/s, 
even though the current may be supplied largely by fewer faster moving particles. 
Zmuda and Armstrong (1974) have reviewed the observations of energetic particles 
which may contribute to the current. Earlier Cole (1963) showed that field aligned 
currents might be studied from the ground using accurate parallactic imaging of 
auroral rays. Such currents cause a twist in the geomagnetic field and a conse- 
quent change in the direction of the geomagnetic field in spuce. It is suggested 
that the parallactic method, calibrated by the backscatter technique of measur- 
ing current, may form a useful and less expensive additional system for ground- 
based study of the currents on a wiser geographical scale. 
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Tho Cufi^ 


In tho vicinity of the polar cusp ionosphere, field-aligned ourront is highly 
structured and a current of about 10” 4 A/m 2 (10" 5 emu) has boen interpreted 
there during a large fluctuation observed by a rocket borne magnetometer 
(I>edley and Farthing 1974). Though the current is probably carried by energetio 
electrons from the solar wind the current doeB imply a mean electron velocity 
(assuming a density of order 5 X 10 1 1 m" 3 at the rocket altitude) of about 
1.2 x 10 m/s along the geomagnetic field 

Asymmetric Dynamo Action 

Od account of asymmetries of magnetic field, wind fields or ionospheric 
electrical conductivity current will flow parallel to the geomagnetic field to 
maintain continuity (Dougherty 1963, Van Sabben 1966, Mishin et al. (1968), 

Maeda and Murata 1965, 1968, Cocks and Price 1969, Stenning 1969, Cole 1971), 
Currents of strengths in the range 5 x 10 G A/m 2 (5 x 10" 11 emu) to 5 X 10” S A/m 2 
(G x 10” l ° emu) have been estimated in these studies. Combined with electron 
densities of order 10 10 m~ 3 (at 1000 altitude) or about 5 x io u m" 3 in the F 
region this implies electron velocities in the range 6 x 10 1 m/sec to 3 x io 4 m/sec. 

These currents should be expected to flow virtually at all times during the 
day on account of the permanent asymmetry of conjugate ionospheric conduc- 
tivities and dynamo fields. 
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Magnetosphere Currents 


By continuity, the field aligned current at the top of the Ionosphere (say In 
the F region) can be related to that In the magnetosphere by 

= h 

where A = cross section of tube of force with unit cross seotlon at the ionospher e 
(A) « Bj /B, Where B denotes magnetic induction and subscript I denotes the 
ionosphere. It follows that along a magnetospheric field line 

— B. = j . = constant 
B 1 1 

or 

v r — rr- 
neBj 

It follows that only in a magnetotube with naB would v be constant. If n falls 
more rapidly with altitude than this, then some advantage in terms of Aw /&j c 
( see equation 7) is gained by studying the currents at higher altitudes. In any 
case such a study should be done in order to try to locate the sourc.*;' region of 
these currents. 

The Pre-Dawn Ionosphere 

Sometimes the ionosphere is sunlit at the conjugate place before sunrise 
on the local ionosphere. This causes heating of the local ionosphere by the 
energy of photoelectrons absorbed in the geomagnetic field tube of plasma. The 
accompanying predawn enhancement in airgiow may be due to thermal excitation 
by F region electrons (Cole 1965) direct excitation (Carlson 1966) or recombination 
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due to electrodynamically induced lowering or downward diffusion of the F region 
(Carlson and Walker 1972) or a combination of those. In the case of lowering of 
the F region by an electric field, field aligned electric current may perhaps be 
involved. Indeed one would expect field aligned current to exist at these times 
on account of the asymmetry of the ionospheric conductivity and dynamo fields 
in the conjugate ionospheres. The measurement of electric current at these 
times should shed further light on these interrelated geophysical phenomena. 

Eclipses 

Eclipses of the sun either locally or at the conjugate point provide oppor- 
tunities to study the field aligned currents associated with the electrodynamic 
asymmetry induced in the ionosphere. 

Within die Ionosphere 

Within the ionosphere winds of different direction at different altitude have 
•different dynamo effects and compensating electric currents parallel to B should 
flow in order to maintain div j_= 0. The study of the scale sizeB and magnitudes 
of these currents from the ground will add immensely to our knowledge of iono- 
spheric electrodynamics. 

Magnetic Substorms 

The interchange of tubes of plasma in the geomagnetic field must be accom- 
panied by field-aligned currents (Cole 1971). This occurs during magnetic sub- 
storms. It will be of interest to study such currents occurring equatorwards of 
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the auroral electrojet where they may be weak and not yet reached by the satel- 
lite magnetometer, because of problems of baseline determination. 

Break-up of Intertropioal red arc 

At night the intertropioal red arc frequently "breaks -up n , causing airglow 
enhancements of irregular shapes and sizes (Barbier 1958, Petersen et al., 1966). 
These are likely caused indirectly by irregular electric fields which cause the 
interchange of tubes of plasma at low magnetic latitudes. These interchanges 
will be accompanied by field-aligned currents (Cole 1971). The arcs occur at 
about ±15° from the magnetic equator. The incoherent backscatter facility at 
Arecibo, Puerto Rico is well placed for the study of the field aligned currents 
that must accompany the break up of the inter tropical red arc. 

Artificial stimulation of the plasma line 

By the artificial introduction of electron beams from a satellite, or space 
laboratory into the ionosphere the plasma line may perhaps be stimulated at 
night at mid-latitudes for subsequent observation by ground based incoherent 
backscatter iacilities. 

Space Shuttle or Orbiting Incoherent Scatter Systems 

It would add greatly to our knowledge of field aligned currents to operate a 
mobile incoherent facility in space, making the measurements suggested in the 
foregoing. 
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CONCLUSION 


It has been demonstrated that the incoherent backsoatter technique may be 
used to study virtually the whole range of geophysical phenomena involving field- 
aligned electric currents. In that it is able to give absolute measures of currents 
stronger than about 10" 7 A/m 2 (10" 12 emu) it has the potential for the greatest 
dynamic range of the techniques presently available to study these currents. 

The existing facilities at Arecibo (Puerto Rico), Saint Santin (France), Millstone 
Hill (USA) and Chatanika (USA) are well placed for many of the studies suggested 
in the foregoing report. The proposed EISCAT facility in Scandanavia as planned 
will be well placed to study the field-aligned currents associated with auroral 
current systems. Improvements in the signal to noise ratio on account of prox- 
imity to the currents would enable the use of lower power systems on a Space 
Shuttle or orbiting satellite. 
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TABLE I 

Statistical Errors Av in m/s (or AJ In Am' 1 for ne ■ 5 x lo 11 m~ 3 
for given Af (see fig. 1) and signal to noise ratios (S/N) and for 
an operating wavelength \ ■ .3m and an integration time t ■ 1000 s. 


25 kHz 

50 kHz 

100 kHz 

Av * 106 m/s 

Aj = .85 x 10"* Am -2 

150 m/s 

1.2 x lCT* Am' 2 

210 m/x 

1.68 x 10" 5 Am" 2 

Av = 10.6 m/s 

AJ = .85 x 10' 6 Am" 2 

15 m/s 

1.2 x 10" # Am' 2 

21 m/s 

1.68 x 10" 6 Am* 2 

A v = 1.06 m/s 

AJ = .85 x 10" 7 Am" 2 

1.5 m/s 

1.2 x l(T 7 Am' 2 

2.1 m/s 

1.68 x 10" 7 Am" 2 

















FIGURE CAPTIONS 


Figure 1. Schematic representation of the plasma lines. 
Figure 2. Diagram showing the inter-relationships of A co/a , 
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Figure 1. Schematic representation of the plasma lines 
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